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“ NASA f s Advanced Exploration Systems (AES) program is pioneering new approaches for 
rapidly developing prototype systems , demonstrating key capabilities , and validating operational 
concepts for future human missions beyond Earth orbit” (NASA 2012). These forays beyond 
the confines of earth’s gravity will place unprecedented demands on launch systems. They 
must not only blast out of earth’s gravity well as during the Apollo moon missions, but also 
launch the supplies needed to sustain a crew over longer periods for exploration missions 
beyond earth’s moon. Thus all spacecraft systems, including those for the separation of 
metabolic carbon dioxide and water from a crewed vehicle, must be minimized with respect 
to mass, power, and volume. Emphasis is also placed on system robustness both to minimize 
replacement parts and ensure crew safety when a quick return to earth is not possible. Cur- 
rent efforts are focused on improving the current state-of-the-art systems utilizing fixed beds 
of sorbent pellets by evaluating structured sorbents, seeking more robust pelletized sorbents, 
and examining alternate bed configurations to improve system efficiency and reliability. 
These development efforts combine testing of sub-scale systems and multi-physics computer 
simulations to evaluate candidate approaches, select the best performing options, and opti- 
mize the configuration of the selected approach, which is then implemented in a full-scale 
integrated atmosphere revitalization test. This paper describes the development of atmos- 
phere revitalization models and simulations. A companion paper discusses the hardware de- 
sign and sorbent screening and characterization effort in support of the Atmosphere Revi- 
talization Recovery and Environmental Monitoring (ARREM) project within the AES pro- 
gram. 


. Nomenclature 

a saturation capacity in Toth equation, mol kg -1 kPa -1 

ao Toth equation parameter, mol kg -1 kPa 1 

a/ superficial free flow area, m 2 

a w column cross-sectional area, m 2 

b equilibrium constant in Toth equation, kPa 1 

bo Toth equation parameter, kPa' 1 

c concentration, mol nr 3 : also parameter in Toth equation, K 

Co inlet concentration, mol nr 3 

c pg gas heat capacity, J kg -1 K 1 

c ps sorbent heat capacity, J kg -1 K 1 

Cpw column wall heat capacity, J kg -1 K 1 

h sg sorbent to gas heat transfer coefficient, W nr 2 K 1 
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axial dispersion coefficient, m 2 s 4 
Toth equation parameter, K 4 

column wall to gas heat transfer coefficient, W nr 2 K 4 

column wall to ambient heat transfer coefficient, W nr 2 K 4 

gas conduction, W nr 1 K 4 

mass transfer coefficient, s 4 

sorbent conduction, W m 4 K 4 

column wall conduction, W nr 1 K 4 

bed height, m 

sorbent loading, mol kg -1 

partial pressure in Toth equation, kPa 

particle Peclet number 

column inner perimeter, m 2 

column outer perimeter, m 2 

average adsorbed concentration, mol m -3 

equilibrium adsorption concentration, mol nr 3 
pellet radius, m 

time, seconds; also heterogeneity parameter in Toth equation 

Toth equation parameter 

gas temperature, K 

sorbent temperature, K 

column wall temperature, K 

inlet temperature, K 

axial coordinate, m 

differential heat of adsorption, kJ mol 4 
void fraction 

void fraction at the column wall 
interstitial velocity, m s 4 
gas density, kg m -3 
sorbent density, kg nr 3 
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p w column wall density, kg nr 3 

D w Distance to nearest wall, m 

G Inverse of distance to nearest wall = 1/D W , nr 1 

g w diffusion coefficient in Eikonal (distance to wall) equation 

pmo Partial pressure of water vapor, Pa 

k Permeability, m 2 

Q Heating source due to sorption, W nr 3 

/3f Forchheimer coefficient, kg m -4 

Ri Reaction rate, mol nr 3 s _1 

C Void scaling 

s Porosity 

Sao ideal perfectly packed porosity (e.g., without the A1 mesh) 

m Bed mass, g 

R Bed internal diameter, mm 

/ Column wall thickness, mm 

Al Langmuir surface area, m 2 g _1 

p Dynamic viscosity, Pa-s 

£m Volume fraction taken up by Al mesh 

rj Sorption efficiency 

II. Introduction 

Predictive simulation tools are being developed to reduce the hardware testing requirements of the Atmosphere 
Revitalization Recovery and Environmental Monitoring (ARREM) project. Although sub-scale testing is required to 
establish the predictive capability of the simulation, the much greater cost of extensive full-scale testing can be lim- 
ited to that required for the confirmation of analytical design optimization studies. Non-recurring costs of predictive 
simulation development are non-trivial, however, once predicative capability is established, geometric reconfigura- 
tion of a model is straightforward. A predictive simulation capability provides numerous additional benefits. Under- 
standing of complex processes is greatly increased since process conditions (temperature, pressure, concentrations, 
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etc.) may be examined anywhere in the adsorption column. Weaknesses in a prototype design can be readily identi- 
fied and improvements tested via simulation. Finally, the predictive simulation provides a powerful tool for virtual 
troubleshooting of deployed flight hardware. 

III. Development of Fixed Bed Modeling Approach 

Adsorption in fixed beds of pelletized sorbents is the primary means of gas separation for atmosphere revitaliza- 
tion systems. For the bulk separation of CO 2 and H 2 O, temperature changes due to the heat of adsorption are signifi- 
cant, requiring the modeling and simulation of the heat balance equations. For columns with small tube diameter to 
pellet diameter ratios, as encountered in internally heated columns, flow channeling along the column wall can have 
a strong influence on overall performance. In some cases, the influence is great enough to necessitate the use of 2-D 
simulations. 

Breakthrough tests, where a regenerated column is challenged with a constant inlet of sorbate and carrier gas, are 
used to determine mass transfer coefficients via empirical correlation. The mass transfer coefficient for a sorbate/ 
sorbent pair may then be used to simulate cyclic, regenerative adsorption processes of interest. The development of 
fixed bed models and verification using breakthrough test results is described below. 

A. Experimental 

Data from fixed bed testing (Mohamadinejad, Knox and Smith 2000) conducted at NASA’s Marshall Space 
Flight Center (MSFC) was used to develop and verify fixed bed models in the COMSOL multi-physics simulation 
package. This test stand had a unique capability where gas concentration was captured not only far downstream of 
the fixed bed but also along the fixed bed centerline in three locations (inlet, midpoint, and endpoint). Temperatures 
were also captured at the same locations as shown in Fig 1. The difference between the exit concentration at the ax- 
ial center the concentration after mixing provides valuable insight into non-ideal flow behaviors such as channeling 
along the wall and axial dispersion, which we will see are difficult to capture in 1-D adsorption models for certain 
sorbate/sorbent systems. Table 1 provides the properties of the adsorbent and fixed bed used to develop and validate 
the appropriate modeling approaches to support AES ARREM CO 2 removal design studies. 



Figure 1. Breakthrough Test Apparatus. The fixed bed test 
section is the center 10 inches of the column. 
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Table 1. Properties of the adsorbent and fixed bed 


Adsorbent 


R p = 1.02 mm 


Fixed-bed 


Pellet radius 
Particle density 
Skeletal density 
Heat capacity 
Langmuir surface area 


p s = 1180 kg nr 3 
psk 2040 kg nr 3 
c ps = 920 J kg -1 K _1 
Al = 463 m 2 g _1 


Bed height 
Bed mass 

Bed internal diameter 
Column wall thickness 
Wall heat capacity 
Wall density 


L = 0.254 m 
m = 396 g 
Ri = 47.6 mm 
1= 1.59 mm 
c pw = 475 J kg' 1 K- 1 
p w - 7833 kg nr 3 


B. 1-D Model 

One-dimensional models are favored for speed and simplicity in fixed-bed simulations. However, in some cases, 
higher dimensions must be used to attain a sufficient accuracy in simulation results. This is illustrated in the results 
sections below. 

1. 1-D Mathematical Model 

Equations implemented in the COMSOL PDE General Form physics mode for the 1-D model are summarized 
below. The following assumptions were used in the derivation of these equations: 

• All mechanical dispersion effects are lumped together with molecular diffusion in the axial dispersion term. 

• Plug flow is assumed, i.e., there is no gradient of velocity, concentration, temperature, or porosity in the radial 
direction. 

• Velocity in the axial direction is not compensated for loss of sorbate since the sorbate gas mole fraction is « 1 . 

• Velocity is temperature compensated via the ideal gas law. 
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Column Heat Bal- 
ance 
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2. Results and Discussion: Breakthrough Tests and 1-D Simulations for Carbon Dioxide on Zeolite CaA 

Figure 2 provides the concentration history and temperature history respectively for breakthrough tests and 1-D 
simulation results with carbon dioxide on zeolite CaA. Note that the mass transfer coefficient k m is treated as a fitted 
parameter. The concentration comparison is quite favorable for the Mixed concentration, which is the desired result 
for a 1-D model. The temperature comparison is also quite favorable; although the peak temperature is slightly 
higher, the point in time when the simulated temperature begins to decrease closely matches test data, indicating 
properly timing of the simulated adsorption front. 

A comparison of the experimental Exit and Mixed concentrations in Fig. 2 reflect the influence of channeling at 
the wall of the column, where packing density is lower allowing greater velocity and earlier saturation. Since the 
average flow rate is higher than that at the centerline of the column, the Mixed data (reflecting the average concen- 
tration for the column) breaks through earlier than the Exit data (reflecting the centerline concentration). 

The plug flow 1 -D model cannot capture cross-flow velocity gradients evident in this data. The axial dispersion 
term was varied from the value representing molecular diffusion to the highest value provided by standard correla- 
tions, yet has negligible effect on the simulation results. Nonetheless, for the C02/zeolite CaA system, the simula- 
tion provides a good representation of adsorption physics and is suitable for use in development studies. 




Figure 2. Breakthrough test data and 1-D simulation results for CO 2 on zeolite CaA. Concentration history 
(left) and temperature history (right). Experimental data are shown as symbols. Simulation data at the inlet 
(2%), midpoint (50%), and exit (98%) are shown as lines. Three values for axial dispersion (units are m 2 s' 1 ) 
are compared in these figures, however, their influence is negligible on simulation results. 
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3. Results and Discussion: Breakthrough Tests and 1-D Simulations for Water Vapor on Zeolite CaA 

Figure 3 provides the concentration history and temperature history respectively for breakthrough tests and 1-D 
simulation results with water vapor on zeolite CaA. Once again the mass transfer coefficient k m is treated as a fitted 
parameter. Additionally, the void fraction was varied for better visual comparison of the test and simulated break- 
through curves. For the FhO/zeolite CaA system, is was not possible to simulate the Mixed concentration data as in 
the C02/zeolite CaA system with the 1-D model. Zeolite CaA has a much greater affinity for water then C02 result- 
ing in much longer time before breakthrough. This provides more time for the effect of channeling to propagate 
down the column as evidenced by the increased gap between the Mixed concentration (diamonds) and Exit concen- 
tration (squares). 

Although the comparison of test and simulated concentrations at the column midpoint and exit is favorable, the 
desired result for design studies is the mixed concentration. The comparison of test and simulated temperatures is 
much less favorable: the later drop in temperature in the simulation indicates that the simulation is not capturing the 
timing of the adsorption process well. 

The data presented here indicates that the 1-D simulation of the FbO/zeolite CaA system is non-conservative and 
not sufficiently accurate for design studies. 2-D axisymmetric models under development for this purpose are de- 
scribed next; these include the variations of porosity near the column wall. 



Figure 3. Breakthrough test data and 1-D simulation results for H 2 O on zeolite CaA. Concentration history 
(left) and temperature history (right). Experimental data are shown as symbols. Simulation data at the inlet 
(2%), midpoint (50%), and exit (98%) are shown as lines. Three values for axial dispersion (units are m 2 s' 1 ) 
are compared in these figures, however, their influence is negligible on simulation results. 


C. 2-D Axisymmetric Model 

The COMSOL 2-D axisymmetric model was constructed using four physics modes. 

1. 2-D Axisymmetric Mathematical Model 

Equations from the four COMSOL physics modes utilized in this simulation are shown below. Detailed informa- 
tion on the origin of these equations may be found in the COMSOL User Guide (COMSOL 2009). 


• Free and Porous Media Flow, including Darcy and Forchheimer terms (Eq. (10)) 
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• Transport of Diluted Species (Eq. (11)) 

^ + V.N = R. 
dt 
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N - - D l Wc + u .c 

• Distributed ODEs and DAEs (adsorbent mass balance; same as Eq. (3)) 

• Heat Transfer (Eq. (12)) 


(/*',. \,j T + h Q 

(p c X Q =( x - e )p. c p. +£ p s c n 

k EQ = ( ! - £ ) k s +£ k g 

Porosity variation is accounted for in Eq. (13), where y is the distance to the wall 
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2. 2-D Axisymmetric Results and Discussion 

Using Eq. (13) to account for the porosity variation near the wall with N = 4, the CCE/zeolite CaA system was 
simulated, with the results shown in Fig. 4 and 5. The porosity variation arises from the lack of pellet nestling near 
the canister wall. It is evident in Fig. (4) that the inclusion of the porosity effect in the simulation results in a faith- 
fully capturing both the mixed and centerline concentrations. This approach is preferred for the CCE/zeolite CaA 
system when higher accuracy is desired, although at the cost of greatly increased computational time. 

It was shown earlier that the 1 -D approach is does not provide acceptable simulation results for the H20/zeolite 
CaA system. The 2-D axisymmetric approach is also being developed for this system; due to the higher numerical 
stiffness of this system, there are convergence issues which must be resolved. 
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Figure 4. Breakthrough test concentration data and 2-D axisymmetric simulation results for CO 2 on zeo- 
lite CaA. 



Time (seconds) 


Figure 5. Breakthrough test temperature data and 2-D axisymmetric simulation results for CO 2 on zeo- 
lite CaA. 


IV. Isothermal Bulk Desiccant Development 

A water-saving bulk drying stage prior to downstream AR trace contaminant control and CO 2 removal processes 
is under development by NASA. Both membrane and fixed bed approaches are under consideration. The primary 
goal is to continuously remove at least 80% and up to 100% of water from a process air stream. Ideally, the bulk 
dryer will operate without active heating during regeneration to minimize power requirements. This may be accom- 
plished by using the dry CO 2 removal system outlet air to regenerate the bulk dryer via a purge which returns the 
captured water back to the cabin atmosphere. Process optimizations such as passive thermal control are of prime 
interest. Development of a bulk drying stage for a water-save AR process at the NASA’s Marshall Space Flight Cen- 
ter (MSFC) has focused on a silica gel-based isothermal bulk desiccant (IBD) system with emphasis on evaluating 
the performance benefits of thermally linking the adsorbing and desorbing beds. 

Experimental efforts to date have included (1) a 4-cell thermally linked system with silica gel sorbent and 17 liter 
per minute flow rate and (2) a 2-cell non- thermally linked system with silica gel sorbent and 140 liter per minute 
flow rate (Perry, Howard, Knox and Junaedi 2009). A mild pressure swing was applied to increase the performance 
of the 140 liter per minute system. 

In the ARREM program, computer modeling and simulation is being used to compare process improvements via 
thermal linking with the additional complexity required. Two analytical efforts are described below; the first in- 
cludes only the thermal balances and arbitrarily imposes the heat of adsorption in order to provide a comparison of 
hardware approaches for thermal linking. In the second effort, adsorption physics are explicitly modeled in 3D using 
an approach modified from the 2D axisymmetric approach previously described. 

A. 3D Thermal Modeling and Simulation - Thermal Linking Design Trade Study 

A trade study has been performed that looks at alternative designs to enhance the heat transfer between adjacent 
desiccant beds in a moisture removal system. The study originally focused on the use of fins running transversely 
between the beds and looked at the effect on performance of variation in fin spacing and thickness. The study was 
subsequently expanded to include designs incorporating mid-bed heat transfer blocks, metal foam and electron beam 
melting (EBM) lattice designs. This study looked solely at the thermal performance of the various designs and as- 
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sumed the adsorbing/desorbing performance was identical in all cases. The finite difference thermal analysis pro- 
gram, Thermal Desktop by Cullimore and Ring, was used to predict the thermal performance of the different bed 
designs. 

1. Modeling and Simulation 

The thermal model was based on 1/3 of a six bed full scale concept, 
incorporating 2 beds each 5.625” deep (in the flow direction), 1.5" wide 
and 7.5” long. Moisture laden air enters one bed at a rate of 20/3 scfm 
and a temperature of 80°F. Simultaneously, dry air at the same rate and 
temperature enters the adjacent bed but from the opposite direction. 

This flow scenario continues for 15 minutes, at which point the flows 
are switched. The simulation continues in this manner, switching every 
15 minutes, until quasi-steady state conditions are achieved. 

The heat of adsorption and desorption were assumed to be equal and 
opposite at a rate of 243.5 btu/hr per bed. This heating(cooling) was 
applied, in its entirety, to l/9th bed sections for l/9th of the half cycle 
time, progressing axially through the bed, in discrete l/9th steps, in the 
direction of flow. Air and desiccant were assumed to be at the same 
temperature. Heat transfer within the bed is assumed to be via the air- 
flow axially through the bed and at the rate of conduction through air, 
radially. Heat transfer to the external environment was ignored. All 
structure is aluminum 6061-T6. 

2. Results and Discussion: Fin Concepts 

Thermal linking performance of fins of different spacing and thick- 
ness was assessed to determine an optimum arrangement. In all cases 
the fins run transversely between beds as shown in the Thermal Desk- 
top model in Fig. 6. Thermal Desktop elements were used in the fins 
and sidewalls so that the thicknesses could be modified without requir- 
ing the configuration to be redrawn. Fin spacings of both 0.5” and 
0.25” and Fin thicknesses from 1/50” to 1/8” were assessed. Interest- 
ingly, the fin thickness had almost no effect on the temperature difference between beds. Reducing the fin spacing 
(increasing the number of fms) has a more significant, but still unsatisfactory effect on the bed-to-bed AT. 

Figure 7 shows average bed temperatures during a half cycle for both beds with no fms and beds with 1/16” fms 
on 1/2" spacing. Decreasing the spacing to 1/4” more than triples the magnitude of the reduction in AT but still re- 
sults in little performance improvement. It is hypothesized that the fin surface area in contact with the flow is much 
too small to result in significant reductions on the bed-to-bed temperature differences and the fin spacing must be 
greatly reduced to increase this area. Unfortunately, very small fin spacing complicates bed packing and leads to 
flow short-circuiting the bed thru the interstitial gaps between the fin walls and the desiccant beads. 



Figure 6. Fin Assessment Thermal 
Model 



3. Results and Discussion: Metal Foam 

Duocel aluminum foam has been used successfully to enhance thermal linking in desiccant beds and was ana- 
lyzed as part of this study. A 1 0 ppi, 8% density foam was chosen for this assessment because it was previously util- 
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ized in subscale tests at MSFC where it demonstrated significant performance in both CO 2 and moisture reduction 
trials. The downside was shown to be an increase in bed pressure drop caused by the reduced bead size required to 
fully load the bed. 

The foam was modeled as porous blocks occupying the same space as the desiccant/air. Thermal parameters 
utilized in the model were from Duocel literature and are reproduced in Fig. 8 (ERG Aerospace Corporation ). 
Thermal performance was vastly superior compared to fins with an average AT of 16°F compared with 31°F for the 
case of no thermal linking enhancement features. Plots of bed temperatures for a half cycle are shown in Fig. 9. 
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Figure 8. Duocel Aluminum Foam Data Showing 
Performance Point Assessed 
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Figure 9. Effect of Aluminum Foam on Tempera- 
ture Difference Between Beds 


4. Results and Discussion: EBM Matrix 

Advances at MSFC with electron beam melting (EBM) of aluminum provide the opportunity to generate deli- 
cate, complex structures within voids. A design of an internal matrix within the desiccant bed was generated in the 
hopes of retaining the overall thermal performance of the aluminum foam but with lower pressure drop and easier 
loading characteristics. The matrix shown in Figure 10 was the result. It consists of offset layers of 0.05” diameter 
strands arranged in 1/2” flat-to-flat hexagonal shapes, and spaced 1/4” apart. The matrix was first analyzed by itself 
to determine its thermal characteristics. This is overlaid on the Duocel foam charts in Figure 11. As is apparent in 
Figure 11, the EBM matrix does not have the density, conductance, or surface area of the metal foam. This is borne 
out in the plot in Fig. 12 where performance is seen to be little better than the non-thermal linking case. 

A design was also analyzed where the strands were increased to 0.625” diameter and the hex size was reduced to 
3/8” flat-to-flat. The thermal characteristics of this design are also shown in Fig. 11 and the thermal linking per- 
formance in Fig. 12. Although the density and thermal conductance are comparable to the metal foam, the surface 
area is much smaller and the performance is degraded because of this. In order to get the performance of the metal 
foam the surface area of the strands would need to be increased. This would lead to something like a star shaped 
cross-section or other complex geometry that have not been pursued. 
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SURFACE AREA 



Figure 10. EBM Matrix Design 

1.5 ” 7.5 ” x 5.625 ” Matrix; 1/2” flat-to-flat dis- 
tance; 0.05” dia. strands; 1/4” spacing between 
layers 



THERMAL CONDUCTIVITY 



Figure 11. EBM Matrix Performance 
Data Superimposed on ERG Duocel Alu- 
minum Foam 



5. Results and Discussion: Heat Transfer Block 

The final concept assessed uses a little different philosophy. Instead of having the thermal linking and the mois- 
ture reduction/bed drying occurring in the same space, the thermal linking is separated out. To accomplish this the 
beds were split in half and a heat exchanger consisting of a solid aluminum block with multiple through holes was 
inserted in the between the bed halves. This concept is illustrated in Fig. 13. The rationale is to allow the tempera- 
ture increase in the adsorbing side and temperature decrease in the desorbing side to occur up to the midpoint in the 
bed then have the heat exchanger essentially reset conditions back to that at the inlets. 

The first design assessed consisted of a 1" thick block with 539 through holes per bed, each 0.1 M in diameter. 
Additional cases assessing the performance of 1/2” and 1/4” blocks as well as the addition of blocks at the inlets and 
outlets were also analyzed. Finally, a case with 8 heat transfer blocks imbedded in the beds, each 3/32” thick was 
analyzed in an attempt to essentially cap the potential performance to be realized with this concept. Performance of 
this case is shown in Figure 14. 
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Figure 14. Performance of Eight In-Bed Heat 
Exchangers 


Figure 13. Heat Transfer Block Design 


6. Conclusions 

All cases analyzed are displayed in Fig. 15 plotted against assumed masses. The eight bed case has significantly 
better performance than the rest, however, manufacturing and loading this configuration may prove challenging. A 
test fixture emulating the 1/3 scale bed layout that the preceding analyses have been based on has been designed and 
is described in a companion paper. It incorporates three 0.2" thick in-bed heat exchangers plus 0.2" thick heat ex- 
changers at the inlets/outlets. While not yet analyzed, it is hoped that this configuration will match that of the 10 ppi 
aluminum foam. The cut-away solid model of this concept is shown in Fig. 16. 
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Figure 16. Cutaway of proposed Test Fixture 
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B. 3D Thermal Modeling and Simulation with 
Adsorption Physics 

1. Modeling and Simulation 

Initial models of the 4-column IBD with an alu- 
minum (Al) foam lattice were conducted used the 
COMSOL Multiphysics code. For the porous flow 
part of the model, the code solves the low Mach 
number, compressible Brinkman equations modified 
by a Forchheimer drag term (Eq. 10 above). 
Forchheimer drag adds a drag term that is propor- 
tional to the square of the fluid velocity, rather than 
just linear, as is the case in Darcy flow, with a coef- 
ficient given by Pf= 1 . 75 p g * V(s / 1 5 0/k) . Heat trans- 
fer is modeled by solving the heat transfer equations 
separately for the porous medium (the beds) and the 
solid Al holder, with appropriate boundary condi- 
tions (Eq. (11) above). For example, the holder has 
a heat flux boundary on its exterior based on the 
ambient temperature of 298.75K. The concentration 
of H 2 O is modeled by solving a diluted species 
transport equation in the form of Fick’s law with 
convection (Eq. (11) above). Since the partial pres- 
sure of H 2 O is -0.01 atm, while the system is at -1 
atm, the dilute transport assumption is valid. The 



Figure 17. Meshed IBD 4-column model. The red 

and blue regions are paired wet/dry inlets/exits of the 
columns. The size of the IBD bed in the three dimen- 
sions are shown in inches. 



time (hrs) 

Figure 18. Temperature comparison. The ‘in’ data (upper left and bottom right curves) are used as boundary 
conditions in the simulation. The values were taken at the center-line of the left red bed in Fig. 1 7, on the inlet 
and exit surfaces. Note the inlet and exit of cell B are spatially fixed, so that the Inlet ’ is where wet air enters, 
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mean loading, q , is modeled, via the Glueckauf approximation, with a simple ordinary differential equation given by 


Eq. (3) above, k is an empirically derived mass transfer coefficient and the equilibrium loading, < 7 *, is a function of 
P and T based on the Toth isotherms. The loading rate is used to determine the heat source in the porous heat trans- 
fer equation, Q , as well as the reactions in the water transport equation, Ru 


Here, Q = — (l — 8 )— — dH and 
dt 


R 


dq (l - 8 ) 
dt 8 


The carrier gas, N2, is modeled as an ideal gas with variable heat capacity, density, thermal conductivity, and 
viscosity. All other physical parameters (A1 thermal conductivity, diffusion coefficient, etc.) are assumed to be con- 
stant except for the system’s porosity and permeability. The porosity, s , varies with location, approaching unity at 
the walls, and is modeled as per the Tobis and Vortmeter prescription (Eq. (13) above) with N= 4, so that: 


1+ 

f-L-il 

e (-4DJd p ) 


£ 
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where Sao is the nominal perfectly fixed bed porosity without an A1 mesh, 


d p is the diameter of the assumed-spherical sorbent pellets, s m is the volume fraction taken up by the A1 mesh, and 
D w is the local distance to the nearest wall. The latter is calculated with a modified Eikonal equation (Fares and 

Shroder 2002) ; VG ■ VG + <7„G(V • VG) = (1 + 2<7„ )G 4 

where G=\/D w and o w =0.1 is a diffusion constant for numerical stability. In this model, inlets and exits are con- 
sidered walls, since the experiment had fine meshes covering the columns, reducing the sorbent packing density at 
the inlets and exits in a way analogous to a wall. The permeability uses the BKZ model (Barenblatt, Zheltov and 

Kochina 1960) that is a function of porosity: K — d 2 £ 3 /(l5C)(l — £ 2 )). The idealized porosity infinitely far from 


any walls is 37% of the volume, with 8 % of the remaining 63% being taken up by the A1 foam (, Soo =0.45 and s m 
=0.08), so that at the walls the porosity is 92%. Note for this small-scale setup, the minimum porosity is 0.385, 
slightly larger than SooSm- The A1 foam, brazed to the canister walls, is, based on Fig. ( 8 ) above, assumed to have a 
thermal conductivity equal to -46% of normal A1 or 92 W m _1 K _1 . The pressure drop across the columns is primar- 
ily a function of geometry and turns out to be quite sensitive to the porosity and permeability models. The above 
models result in a pressure drop of -170 Pa, compared to the experiment value of -164 Pa, while using Sao =0.465 
gives a -140 Pa pressure drop. 

The wet inlet and dry inlet concentrations are held constant at 0.52 and 10 " 10 mol nr 3 , respectively. The outlet 
pressure is held constant at the measured value of 14.2 psi, while the inlet temperature for the dry and wet flow use 
different functions Tm(t) that are half-cycle averages of the experimental data (these Tin functions can be seen in the 
temperature results shown in Fig. 18). 

The inlet boundary conditions for the porous fluid flow are a volumetric flow of 9.998 F/min divided by the 
mean inlet porosity; the flow needs to be normalized in this fashion since the problem starts at the bed inlet and thus 
the interstitial velocity is appropriate. The flow is initialized with a small axial velocity. On the other hand, the ve- 
locity field used by the heat transport equation is the superstitial velocity, so that u = w s. The heat flux between the 
flowing gas and the A1 canister is 16.8 W m-2 K-l, while between the canister and the air it is 1.1 W m-2 K-l. For 
the dilute species transport and heat transport equations, the exits are set to outflow. A ramp for small concentrations 
is applied to the pressure used to calculate the equilibrium loading, q*. Similarly, the inlet concentration is ramped 
up over time. These ramps are applied solely for numerical stability; the results do not change. 

First, a stationary solution is determined first, in order to get good initial conditions for the solver. Then an un- 
steady solution is solved. Every half cycle, the boundary conditions are switched and the unsteady calculation con- 
tinued. COMSOF uses a parallel sparse direct linear (PARDISO) solver that exploits multithreaded parallelism. 

The model discussed here is preliminary and simplistic, since its goal was to determine the applicability of 
COMSOF to the project. For example, a coarse grid, shown in Fig. 17, was used; the calculations are not expected 
to be temporally or spatially converged. Fig. 18 shows the temperature results at the end of 42 simulated 15-minute 
half-cycles compared to the last two half-cycles of the experiment. In the simulation, the wet exit flow initially un- 
dershoots, while the dry exit flow overshoots. Note the inflowing dry air is -1.5K warmer than the inflowing wet air 
and has a different temporal shape. The cause of the -0.2K mismatch between the wet inflowing air and dry exiting 
air and the -0.4K mismatch between the dry inflowing air and the wet exiting air is unknown, but it may be related 
to the fact that the A1 holder serves as a heat sink to the warmer incoming dry air, while the incoming wet air is 
closer to ambient temperature. The dry flow, as it removes water from the silica gel sorbent, cools by -1.4K. The 
wet flow, as it deposits water on the sorbent, heats up by -1.5K. This slight asymmetry is due to the shorter ‘dis- 


American Institute of Aeronautics and Astronautics 

16 



efficiency 


tance of run’ for the hot flow: the wet flow has the whole column for the water to be adsorbed, while the desorbing 
dry flow, since the bed is not fully saturated, only finds water to pick up closer to the exit. 

Fig. 19 shows the water concentration results from the simulation. It is evident that -30 half-cycles are required 
to reach a quasi-steady state; the experiment was run much longer than this. The concentration at the wet inlet, 



Figure 20. Efficiency, t|, of the simulation compared to ex- 
periment over the last half-cycle. During the \ controlled ’ period 
of each half-cycle, the dew point measuring device is calibrating, 
so the resulting concentration and partial pressures are uncertain. 


A 300.57 



T 298.13 


Figure 21. Temperatures (in K) on the bed surfaces at the end 
of the simulation. Counting from left to right, columns 1 and 3 
have wet air flowing downward and columns 2 and 4 have dry air 
flowing upward. 


when it becomes the dry exit, does not drop to zero in the short duration of a half-cycle. The wet flow concentration 
drops by -0.44 mol m-3 before exiting the column, while the dry flow concentration increases by only -0.31 mol m- 
3. This is in-line with the larger temperature change for the wet flow described above. 
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Figure 19. Water concentrations from the simulation. The values are taken from the centerline. The ‘ mid- 
dle ’point is at the axial center of the bed. 


Converting the change in concentration into an efficiency, r|=l-P H 20 0 ut/P H20 in gives the results shown in Fig. X4. 
While the dew point monitor is calibrating (shown in green in Fig. X4), the results are unreliable. Note the function 
used to invert dew point temperature to partial pressure for Fig. X4 is not the same as is used elsewhere in this pa- 
per; thus, the efficiencies are somewhat lower. The simulation has a lower average efficiency overall. The rate of 
increase in the concentration coming out the exit throughout a single half-cycle is also too large in the simulation, so 
that the efficiency at the end of each half-cycle is actually too large; that is, the red curve in Fig. X3 should be flatter 
and and lower. 

The temperatures on the surfaces of the four columns are shown in Fig. 20. Note that the A1 holder (shown in 
outline in Fig. X5) is very nearly isothermal at -299.25 Km with a slight asymmetry of -0.03 K from the desorbing 
to adsorbing sides. Comparing this to the ambient -298.75 K, this is consistent with the above discussion that heat- 
ing slightly dominates cooling. The left-most column in Fig. 21 has wet air flowing downward that is getting heated 
as the H20 gets adsorbed onto the pellets, while the right-most column has dry air flowing upward that is getting 
cooled as the H20 gets desorbed from the pellets. This highlights the importance of accurately measuring the ambi- 
ent temperature, or the holder temperature itself if it is insulated, in future experiments. 

2. Results and Discussion 

The first out-of-the-box attempt at using COMSOL Multiphysics to model the IBD process has resulted in a fa- 
vorable match to data, at least in temperatures. Concentration is not matched as well, but this could be do to a num- 
ber of reasons, from how the data was acquired (compared to the simulation) to variability in physical parameters 
(e.g. sensitivity to the axial dispersion coefficient, DL, used isotropically in the dilute species transport equation, has 
not been investigated). Further work will focus on modeling a thermally unlinked bed, which uses a plastic holder, 
as well as larger scale, higher flow rate experiments. Barring unforeseen issues with those models, COMSOL 
should be a useful tool to explore design space for ARREM IBD. 
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V. Microlith® Adsorber Model Development 

Precision Combustion, Inc. (PCI), via NASA-sponsored projects, has been developing regenerable Microlith®- 
based adsorber modules for the separation of air constituents such as humidity, C02, and trace contaminants to func- 
tion in either closed or open loop operations. The Microlith® adsorbers use a patented technology developed by PCI 
to coat expanded metal with zeolite sorbent crystals to produce a structured sorbent. An advantage of this configura- 
tion is that can be thermally regenerated by passing a current through the electrically-resistive expanded metal. Pre- 
vious developmental work has focused on testing individual Microlith®-based adsorber modules for residual humid- 
ity removal, trace contaminant removal, and C02 removal. The most recent developmental work evaluated the inte- 
grated operation of the residual humidity removal and C02 removal adsorber modules. The adsorber design concept 
is shown is Fig. 22 (Perry et. al. 2009). 


The full scale water Microlith 
“Jelly Roll” was modeled for 
comparison with exit axial veloc- 
ity measurements. The COM- 
SOL Multiphysics modeling 
package, with porous media 
equations and axisymmetric ge- 
ometry, was used. Inspection of 
model results helped to identify a 
probable bypass flow of the 
sorbent bed occurring at the left 
and right ends of the sorbent 
jelly roll as oriented in Fig. 22. 

Testing verified that the previous 
insulation, alumina paper, was 
not thick enough nor compliant 
enough to prevent flow between 
the edge of the jelly roll and the ceramic plate.. 


plate 


Figure 22. Internal cross-section view of the Microlith-based radial 
flow adsorber design concept. 


A. Modeling and Simulation 

Air properties at 20 C and 1.002 bar were assumed for the input and sorbent volumes. Since flow mapping was 
the primary concern, adsorption physics were neglected for the present. Porosity was assumed to be 0.5. Boundary 
conditions at the inlet was designed as standard volumetric rate and at the outlet as no viscous stress. The inlet 
volumetric flow rate was held constant at 600 L/min. Permeability was varied as 9.8e-10 m 2 , 9.8e-9 m 2 , and 9.8e-8 
m 2 for this study in order to match test measurements of overall pressure drop. The model geometry is further de- 
scribed in Fig. 23. 



Ceramic plate & 
aluminum retention 


Ceramic plate 


Custom-made perforated 
tube insert (removable) 
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< Volumes 1 (green) and 2 (blue) 

Air Properties at 20 C, 1.002 bar 
Reactions not analyzed in this study 

< Volume 2 only 

Porous Matrix Properties 
porosity = 0.5 

permeability = sorb_prmblty* 

Edges > 

Axis of Symmetry at red lines 1 and 3 
Inlet at line 2 

standard P&T mass flow rate = mdot_H20** 
Ext Walls and bands at dark blue lines 
Outlet at line 14 
Pressure, no viscous stress 
pO = 0 Pascals 


* varied in this study 

** held constant in this study at 600 L/min 



Figure 23. H 2 O Microlith Model Geometry 


B. Results and Discussion 

The jelly roll exit velocities for the three sorbent permeabilities modeled are shown in Fig. 24. Note that a flow 
rate of zero ft/min was both simulated and measured below 4.7 inches and above 13.2 inches. Immediately adjacent 
to the insulation wrap, flow rates of approximately 50 ft/min were measured. This correlates well with the simulated 
velocities in Fig. 24 using the highest value of 9.8e-8 m 2 for permeability. Figure 25 provides the 3-D velocity map- 
ping of the jelly roll for the lowest and highest permeability values. 
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Velocity magnitude (ft/min) 



Figure 23. Jelly Roll Exit Velocities 



Figure 25. 3-D Velocity Mapping of the Jelly Roll for the Lowest (left) and Highest (right) Permeability 

Values 
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As discussed in the companion design paper, this analysis provided important insight leading to successful cor- 
rective action for the flow channeling problem. Further uses of this model will be to map flow maldistribution fol- 
lowing corrective action. With the addition of adsorption physics, this model will provide a means for optimization 
of cyclic parameters for this particular hardware, and allow for design optimization studies for new Microlith de- 
signs. 


VI. 


CARE Test Facility Vacuum System Analysis 


A. Introduction 

This section describes the fluid analysis of the Common Atmosphere Revitalization for Exploration (CARE) test 
facility vacuum piping system. The CARE concept is to develop a common adsorption canister for applications 
ranging from portable life support systems to long-term habitats. The objective is to reduce hardware development 
costs and increase reliability via long-term testing and operational experience. The vacuum system will be used for 
desorption of CO 2 and water from the fixed bed in a proof of concept test. The purpose of the analysis is to develop 
models to predict the vacuum conductance of the piping system. The resulting models will be coupled for the 2D 
fixed bed analysis. The detailed COMSOL 3D vacuum model was used as a virtual testbed to verify a simplified ID 
pressure drop model using Mathcad. 


c 



Figure 26. Vacuum Piping System 




> 

Closed off 


2500 

B. Experimental 2000 

The inlet and outlet of the piping system are s 1500 

illustrated in Figure 26. The piping starts as 3” fc 100Q 

diameter pipe and increases to 6” diameter after 500 

the valve. The valve internals were simplified and 0 ^ 

modeled to match the test hardware as closely as 0 

possible.. The volume and surface, where the 
pipes bend and are joined, were simplified to 
avoid meshing errors. The 3D model was created 
using SpaceClaim and imported to COMSOL. 

The physics and boundary conditions for the model were determined by a previous test conducted with the same 
vacuum piping and pump configuration. For this test, a controlled mass flow of nitrogen was injected into a fixed 
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Figure 27. Pump Curve 
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bed open to the vacuum system. Various mass flows were injected and pressure and volumetric flow rates were re- 
corded for each case. The vacuum pump curve created from these data is shown in Fig. 27. 


C. Modeling and Simulation 

7. COMSOL Model 

The vacuum flow regime was determined first to identify the appropriate flow equations. For inlet pressures 
greater than 0.064 torr, as expected for CARE test conditions, the Knudsen number was calculated to be less than 
0.01, indicating that the flow is in the continuum regime in the 3 inch pipe. The largest volumetric flow rate and the 
area at the outlet were used to determine the velocity and calculate the Reynolds number of the piping system. The 
highest Reynolds number was Re = 73 indicating that the flow is in the laminar regime. 

The physics used in COMSOL was Single Phase Flow-Laminar Flow- Weakly Compressible. The dynamic vis- 
cosity was defined as a constant of .000017436 Pa*s. An initial condition of 1 Pa was assumed for the entire volume. 

The boundary condition for all walls was assumed to be the no slip condition except for the inlet and outlet. The 
inlet was given a pressure boundary condition which was varied for each case. The outlet condition was interpolated 
from the pump curve data, which was entered into COMSOL as a static pressure curve. 

2. Mathcad Model 

A simplified ID pressure drop model was created in MathCad and correlated with the COMSOL 3D model. 
Density was calculated using the Ideal Gas Law EOS as a function of pressure. The pipe system was divided into 
straight pipes and bends and the head losses were calculated at each component. The equations used for each region 
were derived based on fully developed Poiseuille flow in round pipes; these are summarized below. 
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D. Results and Discussion 

The 3D COMSOL model and the ID Mathcad model showed very similar results for the same boundary condi- 
tions and assumptions. The ID model used compressible equations where possible to achieve as close a match a 
possible to the compressible COMSOL model predictions. 

The outlet pressure comparison from both models is very favorable at lower pressures, though some divergence 
is observed at higher inlet pressures. Figure 28 shows the outlet pressure comparison. Similar results for mass flow 
are observed, with a very favorable comparison at lower pressures, and some divergence at higher pressures. 

The detailed COMSOL 3D vacuum model uses a small amount of computing resources to solve for the cases 
used in this analysis. However, coupling the 3D vacuum system COMSOL model with the 2D fixed bed COMSOL 
model would require a large amount of computing power since the packed bed analysis is complex with multiple 
physics solved simultaneously. Alternatively, a reasonable correlation using the empirical ID equations as shown 
here may be used for the vacuum system, reducing the required computing resources. 
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Figure 29. Outlet Pressure vs. Inlet Pressure 


Inlet P (Torr) 

Figure 30. Mass Flow vs. Inlet Pressure 
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VII. Conclusions 

The need for atmosphere revitalization systems that are optimized with respect to performance, resources, and is 
necessitated by the aggressive new missions planned by NASA. With NASA budgets remaining flat, innovative ap- 
proaches to new system development are required. This paper presents such an approach for the AES ARREM pro- 
ject, where testing is supplemented with modeling and simulation to reduce costs and optimize hardware designs. In 
this paper, we have discussed the empirical determination of mass transfer coefficients using accurate fixed bed 
models in ID and 2D, the optimization of heat transfer for development of a Isothermal Bulk Desiccant (IBD) and 
the application of the fixed bed model in 3D to simulate a cyclic IBD sub-scale test. Studies of the Microlith® Ad- 
sorber flow pattern have been used to troubleshoot performance problems and to obtain a successful solution to the 
flow maldistribution. Finally, the groundwork has been laid for a Common Air Revitalization for Exploration 
(CARE) test by developing the appropriate, simplified vacuum system equations and verifying them against a de- 
tailed 3D multiphysics simulation. 

The efforts represented here will be continued to support the design of Atmosphere Revitalization systems under 
the ARREM project. These modeling and simulation efforts are expected to provide design guidance, system opti- 
mization, and troubleshooting capabilities for atmosphere revitalization systems being considered for use in future 
exploration vehicles. 

VIII. References 

1 NASA. Human Exploration & Operations (HEO). 2012). 

2 Mohamadinejad, H., Knox, J.C. and Smith, J.E. Experimental and Numerical Investigation of Adsorption/ 
Desorption in Packed Sorption Beds under Ideal and Nonideal Flows. SEPARATION SCIENCE AND TECHNOLOGY , 
2000,35(1), 1-22. 

3 COMSOL. COMSOL Multiphysics®. 2009). 

4 Perry, J., Howard, D.F., Knox, J.C. and Junaedi, C. Engineered Structured Sorbents for the Adsorption of Carbon 
Dioxide and Water Vapor from Manned Spacecraft Atmospheres: Applications and Testing 2008/2009. International 
Conference on Environmental Systems (SAE, Savannah, GA, 2009). 

5 Corporation, E.A. Duocel® Aluminum Foam. 

6 Fares, E. and Shroder, W. A Differential Equation for Approximate Wall Distance. Inti. J. Num. Meth. Fid , 2002, 39, 
743-762. 

7 Barenblatt, B.I., Zheltov, I.P. and Kochina, I.N. Basic concepts in the theory of seepage of homogeneous liquids in 
fissured rocks. J. Appl . Math. Mech., 1960, PMM 24, 1286-1303. 


American Institute of Aeronautics and Astronautics 

24 


